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ABSTRACT 

The colliding galaxy pair NGC 2207/IC 2163, at a distance of ^^39 Mpc, was observed with Chandra, 
and an analysis reveals 22 well resolved, luminous X-ray sources. Most of these are ultraluminous 
X-ray sources (ULXs) with Lx ^ 10'^^ erg s^^. The number of ULXs is comparable with the largest 
numbers of ULXs per unit mass in any galaxy yet reported. In this paper we report on these sources, 
and quantify how their locations correlate with the local star formation rates seen in spatially-resolved 
star formation rate density images that we have constructed using combinations of Galex FUV and 
Spitzer 24-/xm images. We show that the numbers of ULXs are strongly correlated with the local star 
formation rate densities surrounding the sources, but that the luminosities of these sources are not 
strongly correlated with star formation rate density. 

Subject headings: stars: binaries: general — stars: formation — stars: luminosity function, mass 
function — stars: neutron — galaxies: individual (NGC 2207/IC 2163) — galaxies: 
interactions — galaxies: nuclei — galaxies: starburst — galaxies: structure — X- 
rays: binaries — infrared: galaxies 



1. INTRODUCTION 

Galaxy collisions can provide information on (i) galaxy 
dynamics, (ii) triggers of star formation, and (iii) the ori- 
gins of ultraluminous X-ray sources, i.e., off- nuclear X- 
ray point sources whose luminosity exceeds the maximum 
isotropic emission expected from a ~10Mq black hole 
(B H): Lx ^, lO'^^ erg s ~^ (hereafter, "ULXs", see review 
by |Feng fc Soria1[201l] ). When one or both of the collid- 
ing gaIa5aes~EiasaTngh gas content, spectacular bursts 
of star formation may be triggered. Some fraction of the 
more massive stars (> 10 Mq) that happen to be formed 
in binary systems evolve to become high-mass X-ray bi- 
naries (HMXBs) with neutron-star (NS) or stellar-mass 
BH accretors. Substantial numbers of UL Xs may also be 
produced. The Antenna e (Whitmore & Schwcizcr 1995 



carried out 



Zezas & Fabbiano"200 2|, the Ca rtwheel ( |Higdoni|1993 
Gao ct al. 2003; Woltcr et al.|2006|), and Arp 147 ( jllappa^ 
port et al.(2010p grandly illustrate these phenomena. Al- 



though the accepted assumption was that HMXBs would 
be Eddington limited (Lx ^ 10'^^ erg s~^), there is a 
growi ng body of evidence linking ULXs to these sources 
'e^ ~ 



20121 



Grimm et al. 20031 iFabbiano 2005 Mineo et al. 



Compelling simulations of colliding galaxies have been 
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1978;|Gerber et al 



pleton et al. 1996 



see, e. g., 'Lynds & Toomre"1976' 'Toomre 
1992; Mihos & Hcrnquist 1994; Ap 



ligdbn fc Wallin 1997j [Struck et al 



2005[ ), and have produced examples that are remarkably 



similar in appearance to the most tidally disturbed col- 
lisional pairs. 

The ULXs seen in these galaxies are of special interest, 
especially those with Lx ^ 10^° erg s^^. Assuming that 
these objects are BHs accreting at the Eddington limit, 
their luminosity gives a limit on the mass of the accret- 
ing BH of ~50 — 100 Mq . These masses are subst antially 
higher th an those of Galactic BHs (M < 20 



Ozel et 



al.|2010[) and could be relate d with the m e tallicity ot the 



host galaxy (e.g., Belczynski et al. 2010 



20101. The population ot ULXs may be t 



Mapelli et al 
leretore consist 
sources of a different nature, depending on their luminos- 
ity. Sources with Lx ^ 10^° erg s^^ might be powered 
by accreting BHs of mass 3 Mq < M < 100 Mq, whereas 
the nature of the brightest sources (Lx ^ 10 erg s~^) 
still represents an enigma. They might be ^10 — 20 Mq 
BHs emitting at ~5 times the Eddington luminosity, 
BHs with masses ^100 Mq emitting at the Eddington 
limit, or the so-called "intermediate-mass black holes" 
("IMBHs" lO^Mf:, < M < 10^ Mq; see, e.g., IColbert 



& Mushotzky 1999) emitting below the Eddington limit 



or some combination thereof. At the highest end of the 
ULX luminosity function (approaching ~10'''"'^ erg sec~^) 
it becomes increasingly difficult to see how the requisite 
luminosity, even if somewhat beamed, co uld be radiated 
near a stellar-mass black hole (see, e.g., Madhusudhan 
et al.|2008 ). The maximum observed luminosi ty for such 
sources is ~10 erg sec" -1 ( [Farrell et al"]|2009| ). 



The IMBHs are of extreme importance as they are 
thought to be the building bl ocks of super-ni assive BHs 
("SMBHs" M - 10^ Mq; e.g., IVolonteri|2010| . However, 
a conceptual problem with IMBH accretors is their for- 
mat ionjjid_rtieiT_subseg^ captu re of a massive donor 



star. 



Portegies Zwart et al. (2004), among others, have 



proposed that runaway star collisions in newly formed 
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massive star clusters lead to the formation of supermas- 
sive stars (e.g., > 500 Mq) which, in turn, evolve to form 
IMBHs. Theoretical problems with this scenario include 
the highly uncertain evolution of supermassive stars, and 
an implausibly high effic ienc y for producing the requisite 
numbers of IMBHs (see [Kmg^ 2004). 
Previous studies have shown a correlation between the 



and the number of luminous X-ray sources (e.g., 


Grimm 


et al. 120031 IMapelh et al.|20101|Swartz et al.|2011 


IJVlineo 


at al.||2012p. The SFRs may be estimated using any of a 



recombination lines (Ha), forbidden lines {[Oil]), far 
infrared (FIR) continuum, thermal radio emissio n lumi 
nosities individually, or i n combination (see , e.g., 'Kenni 
cutt|[l998 , for a review) 



Leroyetal 

a hnear combination of the Galexb'U 



{ 2008[ ) suggest that 
and Spitzer 24 /im 



bands (see their eqs. DIO and Dll) is particularly good 
in this respect. It is possible to follow Leroy et al. and 
utilize this prescription to produce complete SFR images 
of galaxies with a few arcsec resolution. 

In this paper we introduce a new technique to inves- 
tigate the relation between the number and luminosity 
of bright X-ray point sources and the local surrounding 
SFR. This allows us to probe these relations even in the 
case of small numbers of luminous X-ray sources. The 
analysis involves a direct quantitative comparison of the 
spatial structures in these SFR images with the Chan- 
dra X-ray images. In addition, as has been discussed 
(Calzetti et al. 2007; Rappaport et al. 2012, in prepara- 
tion; Mark Krumholz, private communication 2012), the 
Galex intensities tend to indicate the somewhat older re- 
gions of star formation (> 30 Myr, after the obscuring 
dust has already been cleared), while the Spitzer 24 /im 
images reveal younger star formation (i.e., 5 — 10 Myr, 
and still dust enshrouded) which may be more closely 
related to the upper end of the ULX luminosities (i.e., 
with L x > 10^° erg s e c~^). (See also the closely r elated 



I 

h 



work of Sw artz et al. (^20091; Yukita et al. (20101; Kauf 



man et al. 1 2012 )). Ihus, some ot the theoretical ideas 
concerning the formation and evolution of very massive 
binaries can be investigated. 

We report here on our analysis of the relation between 
the location of luminous X-ray sources and the star for- 
mation rate density in the colliding galaxy pair NGC 



2207 & IC 2163 (see, e 
utilize archival Chandra. 



Ga, 



Kaufman et 

Tex. 



aL ^ 

Spitzer aiii 
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rr^o 



We 
Mi- 



cron All Sky Survey (2MASS) images of these galaxies. 
The galaxy pair is at an estimated distance of 39.6^5 *^ 
Mpc (Arnett 1982). This is the redshift-independent 
distance with the smallest uncertainty provided by the 
NASA/IPAC Extragalactic Database (NEDfJ and it is 
based on SN-Ia measurements. A montage of images of 
NGG 2207/IC 2163 taken with HST, Galex, and Spitzer 
is shown in Fig. fl] (NGC 2207 is the larger galaxy on the 
right). Our analysis of the archival Chandra data reveal 
a total of 22 X-ray point sources with Lx ^ 10'^^ erg 
within the D25 ellipse of NGC 2207 & IC 2163. Such a 
production efficiency of luminous X-ray sources per unit 
stellar mass is comparable with that of the Antennae pair 
of colliding galaxies. 



Another advantage of studying the colliding galaxy 
pair NGC 2207/IC2163, in addition to the fact that it 
is relatively close and well studied in numerous wave- 
bands, is that the dynamics o f its collision have been 
extensively modeled (see, e.g., Struck et al.| 2005 and 
references therein). Struck et al. (2005) find that their 
most successful model, which best reproduces the cur- 
rent geometry and morphology of the two galaxies, has 
the following starting conditions. The two galaxies are 
very roughly coplanar as is the grazing collision trajec- 
tory between them. The disk planes of these galaxies 
are inclined only a modest amount (e.g., ~25° — 50°) 
with respect to the plane of the sky. The best model 
results are obtained if the two galaxies had their clos- 
est approach during their first pass when IC 2163 was 
on the western side of NGC 2207 (see Fig.[l| some 300 
Myr in the past. The orbit of the two galaxies is coun- 
terclockwise in Fig.[T] as is the intrinsic rotation of IC 
2163, i.e., the collision is prograde with respect to the 
more compact galaxy. By contrast, the more expansive 
galaxy NGC 2207 is rotating the opposite way, and the 
collision is retrograde with respect to it. The simulations 
show that the spiral arms of NGC 2207 were present be- 
fore the collision, and are not much perturbed by the 
collision. Whereas, the prominent "ocular" feature in IC 
2163 (see especially the Spitzer images) was created by 
the encounter, and its existence apparently shows that 
the collision can not have been going on for more than 
an orbit of the two galaxies. Any star formation that the 
grazing collision has induced has likely been ongoing for 
the past few hundred Myr. The total masses (including 
the dark-matter halos) of the two galaxies were taken to 
be 1.5 X 10" and 1.1 x 10" Mq for NGC 2207 and IC 
2163, respectively. 

2. X-RAY ANALYSIS 
2.1. Data preparation 

We analyzed the publicly available Chandra ACIS-S 
observation of the galaxy pair NGC 2207 - IC 2163, hav- 
ing identification number 11228 and an exposure time 
of 13 ks. The data preparation was done following the 
standard CIACF] threads (CIAO version 4.3, CALDB 
version 4.5.1) for soft (0.5-2 keV), hard (2-8 keV) and 
broad (0.5-8.0 keV) energy bands. The point source de- 
tection was performed in the broad band using CIAO 
wavdetect. We used the v^-series from 1.0 to 8.0 as 
the scale parameter, in order to account for the varia- 
tion of the effective width of the Chandra point spread 
function (PSF) from the inner to the outer parts of the 
analyzed observation, with reference to the Chandra aim 
point. We set the value of the parameter sighthresh 
as the inverse of the total number of pixels in the image 
10~^), in order to avoid false detections. We used 
maxiter = 10, iterstop = 0.00001 and bkgsigthresh 
— 0.0001. The encircled fraction of source energy used 
for source estimation (parameter eenergy) was set to 0.8. 

We computed a monochromatic exposure map for the 
mean photon energy, i.e., 1.25 keV (soft band), 5.0 keV 
(hard band), and 4.25 kcV (broad band), of each band. 
This was done using the CIAO script f luximage, which 
runs the tool mkinstmap to calculate the instrument map. 



^ http://ned.ipac.caltcch.edu/ 



http: / /cxc. harvard. cdu/ciao4.3/index. html 
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Figure 1. Montage of images of NGC 2207/IC 2163 taken with HST, Galex (FUV), and Spitzer (3.6 and 24 /im). The distance to these 
galaxies is ~39 Mpc. 



i.e., the effective-area- weighted exposure map in instru- 
ment coordinates, for the center of each energy band, and 
mkexpmap to calculate the exposure maps in sky coordi- 
nates for each energy band. The exposure maps allow us 
to measure fluxes and hardness ratios for the detected X- 
ray point-like sources, using the photometric procedure 
described below. 

2.2. Source counts 

For each detected point source we measured the count 
rate inside a circular region centered at the source cen- 
tral coordinates produced by wavdetect. The size of 
the circular region for individual sources was determined 
by requiring the encircled PSF energy to be 85% of the 
total PSF energy. In order to do so, the PSF at the 
position of each point source was first constructed and 
then mapped into the World Coordinate System (WCS) 
reference frame of the relative point source image us- 
ing the CIAO tasks mkpsf and reproject_image re- 
spectively. The background regions were defined as an- 
nuli having an inner radius equal to the radius of the 
source region and an outer radius that is 3 times larger. 
The corrected source counts and errors were obtained by 
performing apertu re-corrected photometry , foll owing the 
sa me met hod as in Voss & Gilfanov ( 2007 ) and Mineo et 
aIl([20T2l): 



S 



{b - d)C - dQ 



ha -dp 
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{ha - dpf 



d^a. 



(1) 



(2) 



Here S is the number of net counts from the source, C is 
the number of counts inside the source region and Q is 
the number of counts in the background region, a is the 
integral of the PSF over the source region (expressed as 
a fraction of the total PSF energy), /3 is the correspond- 
ing integral of the PSF over the source and background 
regions, b is the integral of the (effective-area weighted) 
exposure over the source and background regions, and d 
is the exposure integrated over the source region. 

We found four compact sources having background re- 
gions that overlap their neighboring sources. In these 
cases the source count estimation is compromised and 
it was corrected as follows. For the overlapping point 
sources we defined the radius of the circular region that 
included 90% of the encircled PSF. We excluded these 
regions from both the image and exposure map in order 
to respectively subtract the source contribution from the 
background counts and correct the source area. Finally, 
we again performed the aperture-corrected photometry 
described above using the corrected image and exposure 
map. 

2.3. Luminosities and hardness ratios 

X-ray fluxes in the 0.5-8.0 keV band were estimated 
from the count rates measured for each source as de- 
scribed in the previous section. A counts-to-erg conver- 
sion factor was obtained by first extracting the combined 
spectrum of all point sources detected within the D25 
ellipse, except the central source in NGC 2207, which 
may be an active gala ctic nucleus (AGN) partially cov 



ered by de nse clouds (Elmegreen et al. 2006 Kaufman 



et al.|2012 1 . A background spectrum was extracted trom 



large regions located far from detected point sources but 
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NGC 2207/IC 2163 
Chandra 0.5-8 keV 
True Color 





Figure 2. (top) Chandra X-ray image, smoothed with a Gaussian kernel of cr = 2.5" of the galaxy pair NGC2207/IC2163. The red color 
corresponds to the soft (0.2 to 1.5 keV) band, the green color to the medium (1.5 to 2.5 keV) band and the blue color to the hard (2.5 to 
8.0 keV) band. Note that the definition of "soft" and "hard" in the production this figure is different than that used in the data analysis 
and calculation of hardness ratio, (bottom) same as the top panel except that X-ray source numbers from Table 1 and the D25 ellipse are 
shown. 



on the same CCD chip as the galaxy pair under study. 
Some of the background regions were located between 
point sources within the D25 ellipse in order to account 
for the diffuse emission contribution of the galaxy itself. 

Source and background spectra were created using 
the CIAO task dmextract and the associated weighted 
ARF and RMF files were made using, respectively, the 
tasks mkwarf and mkrmf . The average source spectrum 
was binned so as to have a minimum of 15 counts per 
channel and thereby facilitate minimum-x^ fitting. The 



background-subtracted spectrum was modeled as an ab- 
sorbed power law using XSPEC v. 12.7.1b. The best- 
fit model was obtained for uh = 2.8 x 10^^ cm~^ and 
photon power-law index F = 2.1 with = 19-1 for 24 
degrees of freedom. The best-fit power-law index is con- 
sistent with the centroid of the power law photon index 
distribution for luminous X-ray co mpact sources in star - 
forming galaxies, F = 1.97 ± 0.11 ( [Swartz et aL]|2004 ). 
The best-fit column density is a factor of ~3 larger than 
the average Galactic nn in this direction seen in the 
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Leiden/ Argentine/Bo nn (LAB) Survey of Galactic HI 
( [Kalberla et~aI1[2005 1 . 

I'he best-tit spectral model was used to convert the 
count rates of each of the detected point sources into 
fluxes (ergcm~^ X-ray luminosities were calculated 
assuming the distance of 39.6 Mpc ( Arnett,,198"2] ). 

Hardness ratios were used to investigate the spectral 
properties of the sources detected in th e 0.5—8 keV band. 
The procedure described in Sect. 12. 2 1 was applied to the 
reference source list in both soft (d: 0.5-2 keV) and hard 
{H: 2-8 keV) bands. The respective source counts were 
used to calculate the X-ray hardness ratio as: 



HR 



H 



H 



(3) 



The X-ray properties and Chandra positions of the X- 
ray point sources detected above the completeness limit 
are hsted in Table [J 

3. LUMINOUS X-RAY SOURCES 
3.1. Discrete Source Content 

The properties of the discrete luminous X-ray sources 
that our analysis yielded for NGC 2207/ IC 2163 are 
summarized in Tabic [l] In particular, for all sources de- 
tected within the D25 ellipse, we list the source location, 
the net counts after background subtraction in several 
bands, the hardness ratio (defined above), the X-ray lu- 
minosity, and the X-ray fiux. The Chandra X-ray image 
is shown in Fig. [2] 

In all, 28 sources were detected, one of which (source 
no. 7) is li kely a low-luminosity AGN associated with 
NGC 2207 ( [Kaufman et al.||2012[ ), and 6 are just below 
our completeness threshold (as well as below the 'ULX 
limit' of 10'^^ erg s~^; see §3.2). A total of 21 sources are 
sufficiently bright to be ULXs. We note that source no. 
28 has a soft spectrum, and may be part of an elongated 
soft X-ray feature. It is located in the outer spiral arm 
of NGC 2207, ~1.5' N-W from its center at the location 



of the dust y starburst region called feature i (|Elmegreen 
eraL][2000l ). 



In order to show the distribution of the ULX popula- 
tion with respect to the morphological structures of the 
galaxy pair, we plot the X-ray point sources detected 
above the completeness limit, superposed on the HST 
image of NGC2207/IC2163 in Figure [sl The red circles 
indicate the location of individual ULXs, the circle size 
being proportional to the cube root of the 0.5-8 keV lu- 
minosity of the given X-ray source. 

3.2. X-ray luminosity function 



Following the method of Voss & Gilfanov (2006), 
we computed the completeness function K (Lx) ot the 
present Chandra observations of the I?25 region. At 
the assumed distance of 39.6 Mpc, the luminosity, for 
which > 80% of point sources are detected in the 0.5- 
8 keV band, is Lcomp — 1-0 x lO'^^ergs"^ (equivalently, 
K{Lcomp) = 0.8). Twenty two compact sources with 
luminosities above icomp were detected within the D25 
region. Their luminosities range from 1.2 x lO'^^ergs"^ 
to 1.1 X lO'^'ergs"^. After excluding the central AGN, 
we constructed the cumulative X-ray luminosity function 
(XLF); see Figure [4] We modeled the XLF with a single- 
slope power law with a high luminosity cut-off exceeding 



the luminosity of the brightest compact source detected 
in the galaxy pair: 



dN 
dL 



e X 



38 



^38 



L < W^^ erg s~ 



(4) 



A fit of the cumulative XLF using a maximum likelihood 
method yielded a slope of 7 = 2.35lQ'2g. 

This slope is steeper than the slope of ^1.6 that is 
typically found for high-mass X-ray binary (HMXB) lu- 



minosity distributions below ~10 erg s ( Grimm et al. 



20031 
hanc 



3] |Swartz et al.|2011| |i 
a, due to the limitea 



Swartz etal.|2011 IMineo et al. 2012 ). On the other 



due to the limited sensitivity ot this Chandra ob- 



servation we may be sampling only the roll-off of the 
power-law distribution that extends beyond the above 
mentioned slope at lower luminosities to values of Lx in 
the range of 
of 7 



1 — 3 X 10"^^ erg s . A similar XLF slope, 
2.62lg;5Q, is observed in one of the star-forming 



galaxies in the sample of Mi neo et al.| ( |2012[ ), NGC 3079 
for compact sources detected above a completeness limit 
of ~1038ergs-^ 

Interestingly, there is a lack of sources brighter than 
-104°ergs-i in NGC 2207/IC 2163. On the ot her hand, 
ass uming the HMXB luminosity function from |Mineo et| 
al. (2012 1, and rescaling it to match the observed number 
oisources in NGC 2207/IC 2163, the predicted number of 
HMXBs above 10"'" ergs"^ is 2.4 ± 1.5 (assuming a Pois- 
son distribution) , which is fairly consistent with what we 
observe. Luangtip et al. (in preparation) found the same 
evidence of a lack of sources brighter than ^10^° ergs~^ 
in a sample of 17 luminous infrared galaxies (LIRGs) lo- 
cated at distances between ~14 and ~61 Mpc. This may 
be due to the fact that the brightest ULXs should have 
very short lifetimes (e.g., < few xlO^ yr; see Madhusud- 
han et al. 2008), and therefore their numbers and the 
concomitant position of the break in the XLF may be 
dependent on the very recent SFR history of the host 
galaxy. 

We predicted the contribution of background AGNs 
within the D25 ellipse above the completeness lu minos- 



ity based on the work of Georgakakis et aL] ( 2008 ) . Their 
log N — log 5* function for the U.5-10 keV band was con- 
verted to apply in the broad band of 0.5-8 keV. The 
result is that ~1.7 background AGNs with Lx > -/^comp 
are expected to be present and to have a combined lumi- 

l39 arrr CJ - 1 Jjj WOrk 

AGN "luminosities" are computed as 47rD^SAGN(> 



nosityLAGN(> i-comp) = 7.32x10^" ergs" 



-'compy ! 



where D is the distance to the galaxy pair and 
Sagn is the predicted total flux of background AGNs 
above the completeness threshold flux ^comp- The cu- 
mulative luminosity distribution of background AGNs is 
marked in Fig.|4]by a dot-dashed (red) curve. Thus, per- 
haps one or two of the 22 detected sources may actually 
be background AGNs. 

4. STAR FORMATION RATE 

The SFR is one of the most important parameters 
in the investigation of gas-rich galaxies. Over the last 
decade, several studies have demonstrated the existence 
of a tight correlation between the collective number of 
luminous X-ray sources and the integ rated S FR of late- 
type host galaxies (e.g., Grimm et al.[2 003 ; Mapelli et al 

I'his rela- 



2010 Swartz et al. 201 



tion can now be explore* 



Mineo et al.||2012| ). 
in greater detail using spatially 
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Table 1 

NGC2207/IC2163 X-Ray Source Properties 



Source oj2000 <5j2OO0 0.5-8keV 2-8keV 0.5-2keV HR Lx Fx 





(deg) 


(deg) 


(cts) 


(cts) 


(cts) 


(cts) 


/in38 — 1\ 

(10 ergs ) 


/in — 14 —2 - 

(10 erg cm s 


(1) 


(2) 


(3) 


/ /I \ 
(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1 


Q4 08436 


-21, 


.3851 


20 lb 6.1 


5.9 it 4 


14 it 5.4 


— 41 -t 32 


24 it 7.3 


1 26 -t 39 


2 


94 07186 


-21 


.3806 


40 it 8.5 


11 it 5.2 


29 it 7.5 


— CI 46 -t 21 


47 it 10 


2 51 -t 54 


3 


94.07053 


-21 


.3758 


18 ± 6.3 


11 it 5.2 


7 it 4.6 


0.23 it 0.38 


22 it 7.8 


1.18 it 0.41 


4 


94 06623 


-21, 


.3757 


91 ± 12 


30 it 7.5 


61 it 10 


—0 35 -t 1 3 


110 it 14 


5 73 -t 76 


5 


94 07828 


-21 


.3743 


27 ± 7.4 


8.8 it 4.9 


18 it 6.3 


— 35 -t 29 


32 it 8.7 


1 68 -t 46 


6 


94.06943 


-21 


.3743 


20 it 6.1 


2.1 it 3 


18 it 5.8 


—0.80 it 0.27 


25 it 7.5 


1.34 it 0.40 


7t 


94.09183 


-21, 


.3727 


76 it 11 


66 it 10 


9.2 it 5.2 


0.76 it 0.12 


90 it 13 


4.78 it 0.69 


8 


94 07010 


-21 


.3726 


7.3 it 4.4 


3.7 it 3.6 


3.7 it 3.6 


00 -t 69 


9.1 it 5.4 


48 -t 29 


9 


94.11187 


-21 


.3697 


14 it 5.4 


1.2 it 2.8 


13 it 5.3 


—0.83 it 0.37 


21 it 7.9 


1.10 it 0.42 


10 


94.08552 


-21, 


.3696 


27 it 6.8 


6.7 it 4.1 


21 it 6.1 


—0.51 it 0.25 


32 it 8.1 


1.72 it 0.43 


11 


94 07498 


-21 


.3679 


33 it 7.4 


6.6 it 4.1 


27 it 6.8 


—0 60 -t 21 


41 it 9.1 


2 1 9 -t 48 


12 


94 10386 


-21, 


.3641 


7 it 4.3 


2.4 it 3.2 


4.6 it 3.8 


— 32 -t 70 


8.3 it 5.1 


n 44 -t 27 


13 


94.08058 


-21, 


.3642 


14 it 5.2 


2.3 it 3 


11 it 4.9 


—0.67 it 0.39 


16 it 6.1 


0.86 it 0.33 


14 


94.09434 


-21 


.3608 


26 it 6.8 


9.5 it 4.7 


16 it 5.7 


-0.27 it 0.28 


31 it 8.1 


1.64 it 0.43 


15 


94.07195 


-21, 


.3599 


16 it 5.5 


6.5 it 4 


9.9 it 4.6 


-0.20 it 0.37 


19 it 6.6 


1.03 it 0.35 


16 


94.07056 


-21, 


.3527 


29 it 6.9 


7.5 it 4.3 


21 it 6.1 


-0.48 it 0.24 


35 it 8.3 


1.84 it 0.44 


17 


94.10830 


-21 


.3771 


15 it 5.6 


8.1 it 4.5 


7.1 it 4.3 


0.07 it 0.41 


21 it 7.6 


1.10 it 0.40 


18 


94.10480 


-21, 


.3749 


17 it 6.1 


8.1 it 4.6 


9.3 it 4.8 


-0.07 it 0.39 


21 it 7.4 


1.12 it 0.39 


19 


94.07529 


-21, 


.3686 


7.5 it 4.3 


0.94 it 2.6 


6.5 it 4.1 


-0.75 it 0.63 


9.2 it 5.3 


0.49 it 0.28 


20 


94.07520 


-21 


.3648 


9.9 it 5.3 


5.7 it 4.4 


4.1 it 4.1 


0.16 it 0.61 


12 it 6.5 


0.65 it 0.35 


21 


94.11826 


-21 


.3609 


17 it 5.8 


4.6 it 3.7 


13 it 5.2 


-0.47 it 0.35 


24 it 8 


1.28 it 0.42 


22 


94.10091 


-21 


.3865 


8.3 it 4.5 


it 2.2 


8.3 it 4.5 


-1.00 it 0.53 


12 it 6.6 


0.65 it 0.35 


23 


94.08575 


-21 


.3719 


8 it 4.3 


3.4 it 3.4 


4.5 it 3.6 


-0.14 it 0.62 


9.4 it 5.1 


0.50 it 0.27 


24 


94.11424 


-21 


.3712 


6.2 it 4.2 


1.2 it 2.9 


5 it 4 


-0.60 it 0.79 


9 it 6.1 


0.48 it 0.33 


25 


94.11039 


-21 


.3701 


11 it 5.1 


1 it 2.8 


10 it 5 


-0.81 it 0.46 


14 it 6.5 


0.76 it 0.35 


26 


94.10103 


-21 


.3627 


10 it 4.9 


4.5 it 3.8 


5.5 it 4 


-0.10 it 0.55 


12 it 5.8 


0.63 it 0.31 


27 


94.12492 


-21, 


.3790 


30 it 7.3 


7.2 it 4.3 


22 it 6.5 


-0.51 it 0.25 


37 it 9 


1.95 it 0.48 


28t 


94.06592 


-21 


.3673 


6.8 it 4.8 


it 2.3 


7.3 it 4.8 


-1.00 it 0.64 


8 it 5.7 


0.43 it 0.30 



Note. — (1) Source number, (2) Right Ascension (RA), (3) Declination (Dec), (4) - (6) net counts in broad (0.5—8 kcV), hard (2—8 
keV) and soft (0.5—2 kcV) bands respectively, computed with eq. Uncertainties are quoted to 1 o" and were obtained by mean of 

cq. l|2jl. (7) Hardness ratio, computed with eq. Uncertainties were obtained by applying error propagation to the uncertainties in 

the net counts. (8) X-ray luminosity in the 0.5— 8 IceV band, (9) X-ray flux in the 0.5—8 kcV band, f Central Active Galactic Nucleus, % 
Extended soft X-ray source. 




Figure 3. The HST image of NGC2207/IC2163 with the filled red circles superposed indicating the locations of the X-ray point sources 
detected above the completeness limit. The size of each circle is proportional to the cube root of the luminosity of the individual X-ray 
source. 
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Figure 4. Cumulative X-ray luminosity distribution of the com- 
pact sources detected within the D25 ellipse, above the complete- 
ness luminosity of the galaxy pair NGC 2207/IC 2163 (solid, black 
line). The central AGN is not included. The dash-dotted (red) 
curve shows the pre dicted level of resolved b ackground AGNs, 
based on results from |Georgakakis et al.| l |2008[ |. 

resolved images of SFR surface density. Such SFR im- 
ages can be constructed from Spitzer and Galex a rchival 
data, follow ing a recent technique introduced by |Leroy| 
et al. ( 2008 ) . This will allow us to investigate the spatial 
distribution of the X-ray point sources and their lumi- 
nosities as a function of the local SFR at their location. 

4.1. Star formation rate surface density images 



We adopt ed the recipe provided by eq. (Dll) in Leroy 



et al. ( 2008 ) to estimate the spatially- resolved SFR distri- 
EiitIon"m~the NGC2207/IC2163 system. The SFR den- 
sity in units of Mq yr~^ kpc~ 
by: 



was therefore estimated 



SsFR = 8.1 X 10-^ Ifuv + 3.2 x 10-3/24^™ (5) 

where Ifuv a-nd are in units of MJyster^^. The 

combination of FUV and 24 fim images is relatively easy 
to implement because of their reasonably similar angular 
resolutions (4" and 6" FWHM, respectively) and sensi- 
tivities. 

Importantly, this SFR estimator is sensitive to both 
dust-obscured and exposed star formation activity. As 
the FUV emission originates from the photospheres of O 
and B stars, Galex intensities tend to indicate the some- 
what older regions of star formation (rpuv ^ 30 Myr), 
where the obscuring dust has already been cleared. The 
24 fim emission originates from dust grains heated by em- 
bedded young ionizing stars and traces the star forma- 



tio n over timescales, T24^m ~ 5 — 10 Myr (Calzetti et 
[ar]|2007l Rappaport et al. 2012, in preparation; Mark 
Krumholz, private communication 2012). The SFR over 
these shorter timescales may be more closely related to 
the upper end of the ULX luminosities (see Sect. [6]). 



For the first term of eq. ([s]) we used publicly avail- 
able Galex far-ultraviolet (FUV, 1529 A) background- 
subtracted images from the All Sky Surveys (AIS) pro- 
graixj^ These images are calibrated in units of counts 
per pixel per second, and are also corrected for the rela- 
tive instrumental response. The units of the FUV image 
were converted into MJyster"^ prior to combining the 
latter with the 24 fim image. 

The second term in eq. ([s]) was based on a Spitzer 
MIPS 24 ^m Large Field image. We used the "post 
Basic Calibrated Data" products which are calibrated 
in MJvster"^ and suitable for photometric measure- 
mentfp^ We measured the 24 fj,m background in a region 
away from the galaxy, and subtracted it from the image 
before combining the latter with the FUV map. As the 
pixel scales of the 24 /im and FUV images are different 
(2.45"/pix and 1.5"/pix, respectively), we spatially in- 
terpolated the 24 /im image in order to match the better 
resolution of the Galex FUV image. This was done using 
the routine HASTROM, from the NASA IDL Astron- 
omy User's Librar}p] The routine properly interpolates 
without adding significant spatial information (i.e., no 
spatial frequency content beyond the intrinsic resolution 
of the original image). In essence, the 24 /xm image was 
simply oversampled. The resulting SFR surface density 
map is displayed in Figure [5] at the same spatial resolu- 
tion and with the same pixel coordinates as the Galex 
FUV image (but limited by the Spitzer resolution) . 

4.2. Integrated star formation rate 

We calculated the overall SFR of the galaxy pair 
NGC 2207/IC 2163 by integrating eq. ^ inside the D25 
region. The result ing equation corresp onds to the fidu- 
cial method from Leroy et al. ( 2008[ ), in particular to 
their eq. DIO, on which their prescription for obtaining 
star formation rate surface density maps (their eq. Dll), 
is based. 

The FUV and 24 /im terms for the integrated SFR 
were obtained by Leroy et al.| ([2008 ) ass uming a stellar 



as m iCalzetti et al. 



(20071 



initial mass function (IMi< 

(i.e., slope —1.3 for the 0.1 — 0.5 Mq mass range and 
—2.3 for 0.5 — 120 Mq). Using this prescription, we ob- 
tained a total SFR of 11.8 Mq yr~^ integr ated within the 
D25 r egion of the galaxy pair (for the Calzetti et al. 
( |2007[ ) IMF). The corresponding value of Sl-'R based on a 
Salpeter IMF from 0.1 to 100 Mq, would be 18.8 Mq yi'^ 
(the IMF conversion yields a factor of 1.59 difference). 
For comparison, the integrated SFR of the interacting 
galaxy pair NGC 4038/39 (The Antenna e), assuming 
a Sal peter IMF as above, is 5.4M0yr~^ (Mineo et al 
20T2|. 



5. SPATIALLY-RESOLVED A^x - SFR, Lx - SFR RELATIONS 

We investigated the occurrence of the X-ray point 
sources and their luminosities as a function of the lo- 
cal SFR. We started from t he S FR density image ob- 
tained as described in Sect. 14.11 A set of SFR den- 
sity bins with constant logarithmic spacing was defined. 
The source number and their collective luminosity above 

^ http:/ /galex. stsci.edu/GR4/?page=mastform 

-"^ http:/ /irsa. ipac.caltech.edu/applications/Spitzor/Spitzer/ 

■"-^ http:/ /idlastro. gsfc.nasa.gov/ 
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Figure 5. SFR density map for NGC 2207/IC 2163, in units of Mq yr ^ kpc obtained by combining Galex FUV and Syitze r 24 / im 
images according the prescription of Leroy et al. (2008). The small yellow circles mark the locations of the X-ray sources. See Sect. Hllfor 
details. 



icomp (respectively N^{L > Lcomp) and Lx(> ^comp)) 
were assigned to each bin of SFR density according to 
the SFR density value at the position of the source. 
We counted the number of pixels in each bin of SFR 
density and their cumulative area. Knowing the area, 
bas ed on the log N — log S function from Georgakakis et 
al. (2008), we calculated the predicted number of back- 
ground AGNs, Nagn{L > icomp), and their luminos- 
ity, Lagn(> ^comp) above the completeness lumin osity 
threshold, in accord with the procedure used in j |3.2[ 
Typically, this amounts to less than one background 
source per bin. These two quantities were subtracted 
to yield NxiL > Lcomp) and Lx(> Lcomp) respectively 
and the resulting value was divided by the total area in 
kpc^ in each SFR density bin. 

The final values of surface density of X-ray point 
sources (sources/kpc^) and luminosity (erg s~^/kpc^) 
corrected for background AGNs, are plotted against the 
value of the SFR surface density in Fig. |6] Pixels with 
SFR density less than 6 x 10~^ Mq yr~^ kpc~^ were not 
used as they are dominated by background noise. The 
latter was measured in two large regions of the SFR den- 
sity image outside the D25 ellipse. The resulting mean 
values in the two regions are (6.4 ±1.3) x 10^^ and 
(7.3 ± 0.9) X 10"^ Mq yr-i kpc"^ with rms in the range 
of 6.4- 6.8 X IQ-''. 

For comparison with the present measurements, in 
Fig. [6] we also plot the multiple-galaxy-wide average 
A^ x - SFR and Lx - SF R relations for HMXBs obtained 
in 



is equal to our iVx. The dashed line in the top panel of 
Fig. [6] shows their eq. (6), which is slightly non-linear. 

Fig. |6] shows that the global relation between cumu- 
lative number of X-ray point sources and the integrated 
SFR of the host galaxy also holds on local scales (top 
panel). The small-number statistics involved did not al- 
low us to study in any detail the Lx — SFR relation (bot- 
tom panel). A more extensive work by Mineo et al. (in 
preparation) , based on a significantly larger number of X- 
ray sources detected in a sample of nearby grand-design 
spiral galaxies, may yield more useful information on the 
linearity of the latter relation. 

6. SFR DENSITY DISTRIBUTION AROUND ULXS 

We studied the distribution of the SFR densities 
around the detected ULXs, with the ultimate aim of con- 
straining the ULX evolution in relation to the star for- 
mation time scale. The analysis that we discuss below is 
based on the SFR density image shown in Figure [5] 

HMXBs can have large runaway velocities due tokicks 
caused by asymmetric explosions in the formation of the 
compact object. Since it has been su ggested t hat ULXs 
mi ght be bright HMXBs ([ Grimm et al.„2003i [Mineo et 

ocities may also be relevant 



al. 



to 



Mineo et al] ( |2012[ ) (their eqs.(20) and (22)). Here 
A^x = N{Lx ^ 10'^'^ erg s~^) as for th e observed data. 
Since the SFR in the Mineo et al. ( 2012 ) relation is based 
on the Salpeter IMF from U.l to iDO Mq, we first ad- 
justed it to be con sistent with the IMF assu med in the 
Leroy et al. (2008) algorithm (see Sect. 4.2 for details). 
We also c oinpare our results w ith the A^ulx ~ SFR rela- 
tion from Mapelli et al. ( 2010 ), whose definition of A%lx 



2012|), high spatial ve 

LXs. A compact object with a high-mass compan- 
ion can have an average runaway speed of the order of 
~40kms~^, in the case of OB supergia nt X-rav binaries 
or ^15kms ~ ^ for Be/X-ray binaries (Chevalier & Ilo- 
vaisky||1998| ). |Coe| (|2005[ ) measured a maximum average 
speed of ~30 km s~ based on the separation of star clus- 
ters and HMXBs in the Small Magellanic Cloud (SMC). 
Other studies fo und that the limiting averag e speed is 
~50knis~^ (e.g., van den Heuvel et al. 2000, Kaaret et 
2004). The average speed among tliose nested here, 
35kms~^, was adopted as a typical runaway velocity 
of our luminous X-ray sources. Assuming that the latter 
move at this speed over a ~30 Myr lifetime, appropri- 
ate for HMXBs and ULXs, they would be displaced by 



al. 
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Figure 6. Relation between the local SFR density in 
NGC 2207/IC 2163 (Mq yr"! kpc'^) and tlie number density 
(A'x/kpc^, top panel) and luminosity density (Lx/kpc^, lower 
panel) of luminous X-ray sources above our completeness threshold 
in Lx. Both these quantities are corrected for the contribution of 
backgro und AGNs. The st ar formation rate density was computed 
from the |Leroy et al.| l |200 8| algorithm as described in Sect |4.l[ The 
solid curves are the galaxy-wide average A'^x ~ SFR (top pan el) and 
Lx — SFR ( bottom panel) relations for HMXBs obtained in |Mineo| 
|et al.| (|2012l| (their eqs.(20) and (22) respectively). We show tnem 
tor comparison, rescal ed in order to matc h the IMF assumption in 
the SFR recipe from [Leroy et aT] | |2008| . For these expressions, 
Ax and Lx S'^e A, L(Lx > 10'^^ ergs~^), as is the case for the ob- 
servational data. The dashed line in the top panel is the relation 
between the galajcy- wide average numbers of ULXs and SFR by 
[Mapelli et al.| l |2010^ (their eq. (6)), which is slightly non-linear. 
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Figure 7. Distribution of SFR density values surrounding the 
ULXs in NGC 2207/IC 2163. The SFR density was computed as 
described in Sect. K\ in the regions immediately surrounding the 
21 individual X-ray sources in the colliding galaxy pair. From top 
to bottom, the histograms show the distribution of SFR density 
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I I I M ill h 



(rod), around sources with 2x 10'^'' orgs 



< 3x 1039 ergs 



< Lx 

(blue), around sources with Lcomp 

< Lx < 2 X lO^^ergs-i (ma- 
genta) and for all pixels where no sources were detected (black). 
The histogram peaks for the three ULX groups are shifted slightly 
toward higher SFRs as the X-ray luminosity range decreases. The 
trend is driven by the FIR emission. In a more compelling way, 
the regions with no source detections peak at much lower values of 
SFR densities than where the X-ray sources are found. 



'--^Ikpc at end of life with respect to their location at 
birth. At the distance of NGC 2207/IC 2163 this cor- 
responds to a maximum proper motion of ~5.5". The 
latter value was therefore used as the radius of each of 
the fiducial circular regions immediately surrounding the 
21 high-luminosity X-ray sources in which to measure the 
level of local SFR activity. 

We divided the sources into three groups according to 
their X-ray luminosities: icomp < Lx < 2x 10^^^ ergs~^, 
2 X lO^^ergs-i < Lx < 3 x lO^^ergs-i, and Lx > 3 x 
lO^^ergs-^ The three groups include 6, 7 and 8 sources 
respectively. For each group, we created a corresponding 
SFR density image by selecting, from the original SFR 
density map, only the pixels which are covered by the 
fiducial circular regions centered on the sources in that 
luminosity group. In addition, we created a SFR density 
map for the area of the D25 ellipse where no sources were 
detected, by excluding the fiducial circular regions of all 
detected sources. Histograms of the values in these four 
images show the distributions of SFRs in the vicinities 
of the sources in each group and, to act as a control, in 
the regions away from all sources (see Fig. ItI). 

The peaks of the histograms for the three ULX groups 
appear to shift slightly (by at most a factor of ~2) to- 
ward higher SFR densities values as the X-ray luminos- 
ity range decreases. This is opposite to what one might 
have naively guessed. The slight trend is driven by the 
24 /Lim emission rather than the FUV emission, i.e., by the 
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younger (< 10 Myr, Calzetti et al. 2007; Rappaport et al. 
2012, in preparation) star formation tracer of the two in- 
gredients used to construct the SFR density image. How- 
ever, it is difficult to infer the statistical significance of 
the apparent trend because the uncertainties in the SFR 
density estimation are not straightforward to evaluate. A 
vertical dashed line corresponding to the position of the 
peak in the histogram for the brightest sources is plotted 
in Fig. [71 More in line with expectations, the distribu- 
tion for the regions away from the sources peaks at much 
lower SFR densities. This shows that the fraction of pix- 
els having high SFR density, -0.01-0.03 Mq yr'^ kpc"^ 
is much higher in the vicinity of bright X-ray sources, 
than in the field as a whole. This confirms our expecta- 
tion that ULXs tend to be located close to star-forming 
regions. 

Fig. [T] also reinforces the conclusion found in regard to 
the A'^x — SFR and Lx — SFR relations, shown in Fig. [6] 
that while the numbers of ULXs are significantly corre- 
lated with the local SFR densities, the luminosities of 
those sources do not appear to be so correlated. 

7. TOTAL X-RAY LUMINOSITY PER UNIT SFR 

After having subtracted the predicted contribution 
from background AGNs, the collective X-ray luminos- 
ity of the compact source population with Lx > icomp 
within the D25 ellipse is Lxrb(0-5 - 8keV) = (5.33 ± 
0.37) X 10'*°ergs-^ The total SFR integrated within 
the same region (Sect. 4.2) is 11.8 MQyr"^. 

We compare t his re sult with the more extensive results 
of Mineo et al. (2012). First, we estimate the integrated 
SFR of NGC 'imjTC 2163 following the same method 
used in the latter paper, and t hat is based on the recipe of 
Iglesias-Paramo et al. ( 2006| ), which assumes a Salpeter 



IMF from 0.1 to 100 Mp,: SFR = SFR9 



■NUV 



SFR 



IR- 



NUV° is the near-ultraviolet (2312 A) luminosity uncor- 
rected for dust attenuation and IR is the 8 — 1000 /im 



luminosity (see Sect. 5 of Mineo et al.||2012 for details). 
We obtain a total SFR of 23.7 M© yr \ similar to that 
obtained in Sect. |4.2| under the same IMF assumption 
(cf. 18.8 Mq yr~^; SFR estimations usually have uncer- 
tainties of —50%). In general FUV and NUV emissions 
are consistent with each other within 5 — 7%, so the main 
difference between the relation above and the integrated 
SFR measured as described in |4.2| is due to the different 
IR estimator. 

Using the calibration from Mineo et al. ( 2012| , (their 
eq. (22)) we find that the above Sl^'R implies that 
ixRB(> 10^6 erg s^i) = 23.7 • 2.6 x lO^^ = 6.2 x 
10^"^ ergs"^. Taking into account the difference in the 
luminosity limit, by integrating the HMXB luminosity 
distribution from the same paper, we obtain £xrb(> 
-^^comp) = 4.8 X lO"*" ergs~^. The total X-ray luminosity 
measured for NGC 2207/IC 2163 is in agreement with 
the latter result to within ~10%. Analogously, we can 

([20121), 

to predict that Axrb(-^ > L, 



) = i^.y, tor the same 
luminosity threshold. Comparing this with the observed 
number of sources, 19.3, the agreement is rather good 
(the difference between the solid line and the points in 
Fig. [6| is much larger, due to the different method of 
calculati ng the SFR). This supports the suggestive re- 
sults by Grimm et al. ( 2003, ) that the ULX population 



- at least at the bottom end of their luminosity range - 
might be an extension of the HMXB population at high 
luminosities. 

8. DISCUSSION 

In Figures [T][3| and [5| it appears that IC 2163 is forming 
stars more actively than NGC 2207 (note the much more 
extensive dust and gas content in the former), although 
the latter galaxy hosts most of the detected ULXs. A 
comparison between UV and 24 /im emissions (Fig. [IJ 
shows that the latter is more enhanced, suggesting that 
the star formation activity in IC 2163 may be more re- 
cent (< 10 Myr). This might indicate a possible age ef- 
fect in the distribution of ULX s (see e.g., §7.2 and §7.3 in 
Shtykovskiy fc Gilfanov|2005| ) . The specific SFR, i.e., the 
tib K per unit stellar mass, could help us in investigating 
the latter hypothesis. We measured the integrated SFR 
and stellar mass of the two galaxies separately. The area 
of the sky covered by IC 2163 was defined by visually 
inspecting its optical and infrared images. We did not 
use the _D25 ellipse as it includes a large fraction of the 
nearby NGC 2207. The SFR was estimated as described 
in Sect, [jland the stellar mass using the ifg-band image 
from the 2MASS Large Galaxy Atlas (LGAQ assum- 
ing an effective B — V color of 0.67 mag (p rovided by 
HyperLedcT ^ and using the calibration from [Bell & de 
jOTii[([200r 



The integrated SFR of IC 2163 is g.lMg yr-i and 
its inferred stellar mass is 5.2 x 10^° Mq. Five ULXs 
out of 21 are located within its area. This yields a 
Nx/SFR ~ 0.55 ± 0.25 (Mg/yr)"! and specific star for- 
mation rate SFR/M^ ~ 1.75 x lO-^^yr"!. Subtract- 
ing these values from those measured within the en- 
tire £'25 ellipse of the galaxy pair, we estimated for 
NGC 2207 a total SFR of 14.6 Mq yr"! and steUar mass 
1.2xlO"M0. This yields a SFR/M^ ~ 1.22x 10^1° yr"! 
and TVx/SFR - 1.09 ± 0.27 (Mg/yr)"!. The difference 
in specific SFR is not large. The difference in Ax/SFR 
is only marginally significant. The source number in 
IC 2163 is a factor o f 1.4 lower than tha t predicted by 
the calibration from Mineo et al. (2012), while in the 
rest of the merging system it is a factor of 1.5 higher 
than predicted. 

We can compare these properties with those of the 
Antennae (NGC 4038/39). Using the same method as 
above to est i mate the integrated SFR and stellar mass, 
Mineo et al.] ( [2012[ ) obtain SFR = 5.4 Mq yr'^ and M* = 
3.1 X 10^" Mq . Within the same area of sky the same au- 
thors detected 5 ULXs. This corresponds to only some- 
what larger numbers of ULXs per unit mass and per unit 
SFR in comparison with those of NGC 2207/IC 2163. 

The above results for galaxy masses, SFRs, specific 
SFRs, and luminous X-ray sources per SFR are summa- 
rized in Table O 

We tried to identify possible counterparts to the 21 
ULXs hosted by NGC 2207/IC 2163. We cross-matched 
their locations with the Naval O bservatorv Merg ed As- 
trometric Dataset (NOMAD)[^ ( [Zacharias et aLj|2004[ ). 
This catalog is a merger of oata from the HipiDarcos, 

http:/ /irsa. ipac.caItech.edu/applications/2MASS/LGA/ 
http://leda.univ-lyonl.fr/ 
-'^ http:/ /www. usno.navy.mil/USNO/astrometry/optical-IR- 
prod/nomad 
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Table 2 

Summary of Star Formation Rates 



Galaxy 
(1) 


r tot 

(10""' ergs s-i) 
(2) 


M, 
(10"AfQ) 
(3) 


SFR 

{Mq yr-i) 
(4) 


SFR/M. 
(10-10 y,-i) 

(5) 


Afx/SFR 
(6) 


NGC 2207 


4.3 ± 0.3 


1.2 


14.6 


1.22 


1.09 ± 0.27 


IC 2163 


1.0 ± 0.2 


0.52 


9.1 


1.75 


0.55 ± 0.25 


Entire D25 


5.3 ± 0.4 


1.7 


23.7 


1.39 


0.89 


Antennae 


1.7 


0.31 


5.4 


1.74 


0.93 



Note. — (1) Galaxy; (2) collective 0.5—8 kcV luminosity of the compact source population 
with with Lx > 10^^ ergs s~^; (3) mass in stars; (4) total star formation rate, SFR — 
4.6 X lO^^^LiR + 1-2 X 10~'^^LNUv.obB as described in §7; (5) specific star formation rate; (6) 
number of luminous X-ray sources per unit SFR. 



Tyclio-2, UCAC-2 and USNO-Bl catalogs, supplemented 
by photometric information from the Two Micron All 
Sky Survey (2MASS) final release point source catalog. 
We used a match radius of 2" bearing in mind that the 
99% uncertainty circle of Chandra absolute positions for 
sources within 3' of the aimpoint is 0.8". The choice of 
2" match radius is rather conservative considering that 
the typical mean error on coordinates in the NOMAD 
catalog is also < 0.8". We found possible counterparts 
in the USNO Bl.O catalog for 7 ULXs, i.e., source nos. 
2, 4, 7, 14, 15, 25, and 28. For all these objects, the cat- 
alog provides, among other measurements, the B- and 
R-band magnitudes and the proper motions. All the op- 
tical objects are separated from the corresponding Chan- 
dra ULX positions by angular distances ranging from 
1.2" (for source #14) to 2.0" (for source no. 4). Based 
on the accuracy of the absolute Chandra positions, this 
indicates that there are no compelling candidate coun- 
terparts among the nearly matched sources. Moreover, 
sources 2, 4 and 15 have measured proper motions, which 
indicate that these are most probably foreground objects. 
A more accurate search of the ULX counterparts would 
require the use of HST data, which are publicly available 
in the HST archive. A forthcoming paper will include 
such a study. 

The search for optical counterparts can also be useful 
to identify possible background AGNs among our ULX 
sample. In particular, sources 4, 16, 21 and 27 are lo- 
cated outside the spiral structures of the two galaxies. 
For the latter three there is no counterpart in NOMAD. 
As mentioned above, one match was found for source 4, 
with a separation of « 2.0" from its Chandra coordinates. 
This source is the brightest X-ray object detected in the 
vicinity of NGC 2207/IC 2163. The R-band magnitude 
of the optical object is 14.7. This can be converted into 
a flux /r ~ 
turn, log(/o.5 



1.7 X 10- 

-8kev/ /r) 



erg cm ^ s ^ which yields , in 



Bauer et 



2.4. According to ^ 

al] ( |2004| this object should then not be an ALn. i'his 
IS further supported by the fact that the hardness ra- 
tio of source 4 differs from that of the central AGN in 
NGC 2207 (see Table [l|, and also by the fact that the 
predicted number of background sources above the lumi- 
nosity of source 4 is only 0.09. Assuming a Poisson dis- 
tribution, the probability that the brightest ULX in our 
galaxy pair is a background AGN is only 8%. Assuming 
that the observed XLF of the ULXs would extend to the 
highest luminosities with the same slope, the predicted 
number of X-ray sources above the same threshold would 
be 1, thus, a ULX interpretation is the most likely. 



et al. 
density. 



( |2008| , we 
a units of 



9. SUMMARY AND CONCLUSIONS 

We have introduced a new technique to investigate the 
spatial and luminosity distributions of X-ray binaries in 
star forming galaxies as a function of the local SFRs. We 
have applied this technique to study the population of 
21 ULXs in the colhding galaxy pair NGC 2207/IC 2163. 
This is comparable with the largest number of ULXs per 
unit mass in any galaxy that we know of, in particular 
the Antennae. 

Using the prescription by iLeroy 
constructed an image of SFlf 

Mq yr-i kpc-^ of NGC 2207/IC 2163, by combining the 
Galex FUV and Spitzer 24-/im images. We find that the 
global relation between the cumulative number of X-ray 
point sources and the integrated SFR of the host galaxy 
also holds on local scales. We were not able to investigate 
in detail the corresponding local Lx — SFR relation due 
to the small number of X-ray sources. 

We studied the distribution of the SFR density around 
the detected ULXs, with the ultimate aim of constraining 
the ULX evolution in relation to the star formation time 
scale. We find that the peaks in the SFR density dis- 
tributions around ULXs appear to shift slightly toward 
higher SFR density values as the X-ray luminosity range 
decreases. The regions with no source detections, how- 
ever, do peak at much lower SFR densities than where 
the X-ray sources are found. The fraction of pixels having 
high SFR density, ~0.01 - 0.03 Mq yr^^ kpc-^ is higher 
in the vicinity of bright X-ray sources, than in the field. 

We find that the number and luminosity of ULXs per 
unit SFR are in agreement, to wit hin 10%, with thos e 
predicted by the global relatio ns of Mineo et al.l (120121). 
This supports the suggestion of Grimm et al. ( 2003 ) that 
the lower luminosity end of the ULX population (e.g., 
1 — 3 X 10'^^ ers s"^) might be an extension of the HMXB 
population to higher luminosities. 

We attempted to investigate possible age effects 
in the distribution of ULXs across the galaxy pair 
NGC 2207/IC 2163. We find that the difference in spe- 
cific SFR (SFR/M^) between NGC 2207 and IC 2163 is 
not large. The difference in iVx/SFR is only marginally 
significant. The source number in IC 2163 is a fact or of 
1.2 lower than th at predicted by the calibration from|Mi-] 



neo et al. (2012), while in the rest of the merging system 
it is a factor of 1.4 higher than predicted. 

We tried to identify possible counterparts to the 21 
ULXs hosted by NGC 2207/IC 2163. In particular, our 
search of the USNO-Bl catalog emphasized possible op- 
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tical counterparts to the brightest detected source (no. 
4), as weh as for sources 16, 21 and 27, which are located 
outside the spiral structures of the two galaxies. In gen- 
eral, the relatively large separations (i.e., > 1") between 
possibly interesting stellar images and the Chandra co- 
ordinates suggest that there are no compelling candidate 
coimtcrparts. The closest optical object to the brightest 
detected source (no. 4) found in the USNO-Bl catalog 
does not have the properties of a background AGN. A 
more definitive search for ULX counterparts would re- 
quire the use of HST data, which is publicly available in 
the HST archive. We plan to pursue this in a forthcom- 
ing paper. 

The authors are grateful to Mark Krumholz, Adam 
Leroy and Pepi Fabbiano for helpful discussions. We 
made use of Chandra archival data and software pro- 
vided by the Chandra X-ray Center (CXC) in the ap- 
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Propulsion Laboratory, California Institute of Technol- 
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sis Center/California Institute of Technology, funded by 
NASA and the National Science Foundation. Helpful in- 
formation was found in the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, un- 
der contract with the National Aeronautics and Space 
Administration. SR, DP, and AL acknowledge support 
from Chandra Grants G01-12111X and GO2-13105A. 
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